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patent application serial number 10/165,556 filed June 7, 2002. AH of the above 
applications are herein incorporated by reference in their entirety as if fully set forth herein. 

Field OF THE Invention 
This invention pertains to the field of medical imaging, particularly to measuring the 
size and shape of organs fi:om three-dimensional images. 



Background OF THE Invention 
Various kinds of methods have been used or proposed in the literature to measure the 
size, shape, and internal volumes of specific organs using ultrasound segmentation analysis 

10 from three-dimensional ultrasound images. Examples of body structures undergoing 

. . ■ • • ' ■ • . ■ • ' . • ■ . ' 

segmentation analysis include heart chambers, psirticuiarly cardiac chambers (left and right 
ventricles, left and right atriums), the prostate, the bladder, and aninio 

The measurement of he^ chambers is described, for example in U.S patent 6,346,124 

- ^. ' ■ ' ' " - '■ . " • 

to Geiser, et ah (Autonomous Boundary Detection System For Echpcardiographic Images). 

15 Similarly, the measurement of bladder structures are covered in U.S. patent 6,213,949 to 
Ganguly, et al. (System For Estimating Bladder Volume) and U.S. patent to 5,235,985 to 
. McMorrow, et al., (Automatic Bladder Scanning Apparatus). The measurement of fetal head 
structures is described in U.S. patent 5,605,155 to Chalana, et al., (Ultrasound System For 
Automatically Measuring Fetal Head Size). The measurement of fetal weight is described in 

20 U.S. patent 6,375,616 to Sofennan, et al. (Automatic Fetal Weight Determination). 

. Most of these techniques are not ideal in that they cannot be generalized beyond the 
application for which they were desig;ned. They use shape or size assumptions that are not 
necessarily satisfied for different applications. An attempt to make a generalized model is 
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described in U.S. patent 5,588,435 to Weng, et al. (System And Method For Automatic 
Measurement Of Body Structures) but is limited in its application due to method coniplexity 
and operational robustness. \.. 

The techniques commonly employed iii segrnentation use active contour modeling and 
5 live wire methods to delineate an object shape within electronic-based iniages, and partial 
differential equation (PDE) based processing of image signals to reduce background noise 
speckle. These three methods have their advantages and disadvantages. For example, the active 
contour modeling method or "snakes" is a region^based boundary method and requires specific 
shape assumptions and the use of deforrnable models to work properly (M. Kass, A. Witkin, D. 
10 Terzopolous, "Snakes: Active Contour Models,'' International Journal of Computer Vision, pp. 
321-331, 1988). However, "snake" niethods require an initial contour that needs to be very 

close to the final desired boundary. In addition, most of these niethods are iterative in nature 

- \ , - ' ■ . . " ^ • - - - ■ . . 

and labor intensive. . : . : ; . 

Live wire methods, championed by A.X. Falacao et al. in "User-steered Image 

. ' • • ' ' , - ' ■ " ... " • • " . ' 

15 Segmentation Paradigms: Live wire and Live Lane," (A.X. Falacao, J.K. Udupa, S. 

Samarasekara, and S. Sharma, Graphical Models and Image Processing, 60, 233-26-^, 1998) 

and E.W. Dijkstra, in "A note on two problems in connection with graphs," (Numerical Math, 
, vol. 1, pp. 269—271, 1959), requires a user to actively provide a starting point and an ending 

point via a mouse, pointer, or equivalent delineation means. Live wire is an interactive 
20 segmentation method which uses a minimum cost path between the last clicked user point to 

the current mouse clicked location. The live wire methodology of Falacao is labor iritensive 

■ . , , ^ , L ■ » 

- -1- .I 

since it requires an iterative approach and generates a series of mouse delineating end-points. 
Furthemiore, live wire methods are limited in accuracy because the simplest optimal path 





-3- 



25315 



PATENT TRADEMARK OFFICE 



DXUC-l-1021AP.doc 



Black Lowe & Graham 

816 Second Avenue 
Seattle, Washington 98104 " 
206.381.3300 • F: 206.381.3301 



PLLC 



contour connecting any two mouse determined end-points is not sufficient in ultrasound images 
that are prone to mouse pointer backtracking. That is, the "user-steering" methods of Falacao^ 
require that the user backtrack the mouse whenever the drawn contour hne has lost track of the 
object boundary as presented on the ultrasound image. 
5 Most ulti'asound processing methods involve some pre-processing techniques to reduce 

" - . " - " - r 

speckle noise and enhance images, such as median filtering methods and PDE-based image 
"... processing. PDE-based methods have the added advantages of preserving edges while 
enhancing images (see P. Perona and J. Malik, "Scale-space and edge detection using 
aniostropic diffiision," IEEE Trans. Pattern Analysis and Machiiie Intelligence, vol. 12, July 

10 1990, pp. 629 - 639; J. A. Sethian, Level Set Methods and Fast Marching Methods, 2nd 
Edition, Cambridge University Press, 1999; S. Osher and L.I. Rudin, "Feature-oriented image 
enhancement using shock filters," SIAM Journal of Numerical Analysis, vol. 27, pp. 919-^ 
940, August 1990; and U.S. patent 5,644,513 to Rudin, et al. "System Incorporating Feature- 
Oriented Signal Enhancement Using Shock Filters"). These methods, especially shock filters 

15 for image sharpening or blind deconvolution, have been used reducing noise speckle for 
computerized-tomography (CT), magnetic resonance (MR),, and positron emission-tomography 



(PET) images. However, shock filters implemented with unoptimized algorithms have been 
found deficient in that the unoptimized algorithms increase or enhance noise speckle, thereby 
degrading the CT, MR, and PET images. 



20 



In iiitrasound images, chamber-like structures (for example: bladder, heart left 



ventricle, prostate, amniotic sac) are further prone to image degradation due to ultrasound 
echoes reflecting from other surfaces parallel to do organ wall boundary. The echoes from the 
parallel surfaces present as bright regions overlapping with the dark chamber boundary walls. 
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Thus, ultrasound-echoing overlaps creates missing lateral boundaries of the organ walls. 

Thereby, inaking it difficult to accurately determine boundary . wall locations where 

segmentation lines need to be drawn. 

Thus, there is a need for a generalized image analysis method to accurately measure the 
5 size, shape, and volumes of : fluid-filled and non-fluid filled organs from analysis of one- 

dimensional (ID), two-dimensional (2D), and three-dimensional (3D) images fi'pm ultrasound, 

CT, MR, and PET procedures. The generalized image analysis method requires optima;l 
\ algorithrhs to minimize background noise speckle and maximize structure segmentation or 

delineation by minimizing the creation of missed lateral boundaries of organ structures. The 
10 generalized image analysis method needs to be broadly adaptable and simple to implement. 



Summary OF THE Invention 

• ^ ■ * * ■ * " . 

The invention is a method having a plurality of segmentation algorithms and steps 
optimized to robustly locate and measure the volume of a fluid filled or non-fluid filled 

15 structure or organ from primary images obtained from instruments probing with or receiving 
and processing electromagnetic and non-elecfromagnetic radiation including ultrasound, CT, 
MR, PET, and radioisotopes used in nuclear medicine procedures that produce a plurality of 
ID, 2D, or 3D images. The niethod is applied to ID scans to obtain resultant 2D primary 
images. Each 2D primary image is made visually enhanced by the plurality of segmentation 

20 algorithms and steps of the method. The method then includes the creation and placement of 
accurately drawn boundary line contours that more clearly delineate or segment the structure 
walls within the 2D image. Upon placing the boundary line contours, the method calculates the 
. : volumes for each delineated structure or organ presented in each primary image. Then, as an 
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aggregate, the method provides the ability to calculate the 2D area and the 3D volume of the 
Structure or organ, fluid filled or non-fluid filled, from the summation of each 2D image. 

The method employs subjecting each primary image to a pre-filtering step using non- 
linear filters to enhance the image. The pre-filtering step includes an image-smoothing step to 
5 obtain maximum contrast between organ wall boundaries followed by an image-sharpening 
step to establish image enhancements. The pre-filtering step has a first algorithm and a second 
algorithm, each algorithm having, a partial differential equation (PDE) formula in non-linear 
form. The first PDE algorithm serves as a heat filter, and the second PDE algorithm serves as a 
, shock filter. The sequence of image enhancement is application of the heat filter followed by 
10 application of the shock filter. Each image is subjected to a combination of the heat and shock 
filters to reduce each image's signal noise a.nd to sharpen each image signal while at the same 
time preserving signal discontinuities. Each. primary image subjected to the heat and shock 
prefiltering algorithms is formed into a second image having a cleiarer demarcation along the 
boundary of the structure or organ walls. 



15 



Each second image is then processed to determine where to locate initial edge points by 



either of two locating miethods. A first locating method is a manual, user-deployed radial point 
method. The radial point method is comprised of lines that are drawn in substantially 30-degree 
increments for establishing where maximum spatial gradients exi$t during transitioning 
between light and dark areas. A second locating method to set initial boundary edge points is 
20 an automatic method to determine the leading and trailing regions of wall locations along an A- 
mode first dimension (ID) scan line that is regulated by a hysteresis threshold value. 

An optimal pathway to connect between the initial edge points is then determined using 
a third algorithm to determine aii. optimal pathway. The optimal pathway between any two 
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points is one defined to have the least cost, where the least cost represents the sum of three sub- 
equations; namely, the edge distance cost, the path direction cost, and the previous contour 

distance cost. The optimal pathway represents an automatically implemented way to draw 

■ ■> ' ' • ■ • " ' . . - . ■ . . ■ 

contour boundary lines along wall structure using the automatically selected initial edge points 

5 from the second locating method or as maritially selected by the radial method. The advantages 

' » . 'i ' ' ' • ■ 

compared to the manual live wire methodology is fewer initial points are needed to establish 
boundary contours, the minimum cost-optimal pathway algorithm is more robust, and does not 
require manuaHnteractive mouse pointer setting and reselecting; 

After the boundaries of each structure or organ are ; determin^^ 
10 thQ 2D area and 3D volume of the organ or structure is determined. For fluid-filled organs, this 
includes but is not limited to biadders, amniotic sacs, prostates, cardiac chambers (left and right 
ventricles, left and right atriums), ovarian cysts^ arteries, and veins. Similarly, non-fluid 

containing organs or structures includes but is not limited to kidneys, kidney stones, gall 

• • ■ ■* ' - . ■ > ■• • " . • . 

bladders, gall bladder stones, and ovaries. / 



15 



Finally, the method may be done remotely using a web server to analyze scan data sent 



from ultrasoundj GT, MR^ and PET devices over intranets or internets vyhere the device 
architecture permits image processing on the server-side of the invention. 



Brief Description of the Drawings 



20 



Figure 1 is an overview of the segmentation algorithm method; 



Figure 2 is an algorithm processing oh a real bladder image presented in a five-panel 



senes; 



Figure 3A is a depiction of a ID scans line forming a 2D scan plane; 
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Figure 3B is a depiction of a plurality of scan planes forming a 3D image cone; 

Figure 4 is a three panel series of 2D images smoothed with varying heat filter 



algorithm iterations; 



Figure 5 A is a graph showing the effects of a shock-filtering algorithm of a blurred 



5 signal; 



Figure 5B is a graph showing the effects of heat and shock filtering algorithms of a 



blurred noisy signal; ■ 

Figure 6 is a pixel-processing algorithm for finding the direction cost; 
Figure 7 is a 24 panel series of segmentation algorithms applied to 24 image planes of a 
10. bladder phantom, each image plane showing the resultant boundary contours; 

. Figure 8 is a 12 panel series of segmentation algorithms applied to a first human 
bladder scanned with 12 images planes, each image plane showing the resultant boundary 
contours;and . V 

Figure 9 is a 12 panel series of segmentation algorithms appUed to a second human 
15 bladder scanned with 12 image planes, each image plane showing the resultant boundary 



contours. 



Detailed Description of the Preferred Embodiment 



Figure 1 shows an overview of a segmentation algorithm 10. The segmentation 
algorithm 10 begins with the receiving of an input image 20, to which a heat and shock 
20 filtering algorithm 30 is applied to the input image 20. Thereafter, a spatial gradient 
algorithm 34 is applied to the image that has been heat and shock filtered. At this juncture, 
initial edge points are located in a point7positioning algorithm 38 to the boundary of the 
organ structure within the image. Thereafter an optimal path algorithm 42 is applied to 
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determine the best path to connect the initial edge points assigned by the point-positioning 
algorithm 38. This optimal path algorithm 42 also uses the result of the spatial gradient 
algorithm 34 as an input. For images already having a drawn contour, a previous image 
contour input 46 is inputted to the optimal path algorithni 42. TThereafter, a compute volume 

5 algorithm 48 is applied to the image, and thie computied volume is reported as an output 
volume52. . *; 

Figure 2 shows the segmentation algorithm applied to a sample image. An original 
bladder image 204 acquires more contrast as shown in a bladder image 206 once subjected to 
the heat and shock filtering algorithms. Thereafter, spatial gradients are depicted as shown in 

■10 a bladder image 208. Initial points are placed as shown in a bladder image 210 when the 
pointrpositioning algorithm is applied. The point-positioning algorithm includes the radial 
method. Finally, a contour along the boundary where the points have been connected by 
applying the optimal path algorithm is shown in a bladder image 212. 

Figure 3A depicts the ID and 2D operational aspects of a single element transducer 

15 that rotates mechanically about an. angle O. A scanline 314 of length r migrates between a 
first limiting position 318 and a second limiting position 322 as determined by the value of 
the angle creating a fan-like 2D scan plane 310. In the image sho\yn in Figure 2, the x- 
axis is the tilt angle ^, and the y-axis is the scanline depth r. The single element transducer 
operates substantially at 3.7 MHz frequency and creates an approximately 18 cm deep scan 

20 line 314 and migrates within the angle O having an angle value substantially near 0.07 
radians. A first motor tilts the transducer 60° clockwise and then counterclockwise forming a 
. 120° sector image. / ; / 
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The 3P ultrasound data is depicted in Figure 3B that shows the scanplane 310 among 

. ' • • ■ ■ ' '. . " ■ ' — . ' ■ " . " . - . ' ■ 

a plurality of scanplanes substantially assembled into a rotational array. After gathering a 2D 
sector image, a second motor rotates the transducer 15° to gather the next 120° sector image. 
This process is repeated until the transducer is rotated through 180°, resulting in a cone- 
shaped 3D image 340 data set with 12 planes. The 3D image 340 is a rotational array of 
assembled scanplanes. The tilt angle indicates how much the scanline is tilted from the 
center line in 2D sector image, and the rotation angle 0, identifies which rotation plane the 
sector image lies in. Therefore, any point in this 3D. data set can be isolated using plane 
coordinates expressed as three parameters, \P(r,^,^). The rotational scanplane array as 
depicted by 3D image 340 in FIGURE 3B represents one preferred embodiment. Other 



preferred embodiment includes scanplanes arrays assenibled from wedge scans and 
substantially rectangular translational scans. Yet other preferred embodiments include images 
formed by line arrays, either spiral distributed or reconstructed random-lines. The line arrays 
are defined using points identified by the coordinates expressed by the three parameters, 

* • ■ m ' ~ • ' * ' 

P(r,^,^), where the values or r , ^ , and (9 can vary. 

The effects of applying the heat filtering algorithm is shown in Figure 4. The effect of 
10 iterations of heat algorithm is shown in image 410. The effect of 50 iterations is shown in 

image 420. The effect of 100 iterations is shown in image 430. With more iterations of heat 

■ •' ' - ■ ■ . . " ' ' . . 

filtering algorithms, the more the iniage beconies blurry. 

For scanplanes substantially equivalent to scanplane 310, the algorithms are expressed 

in two-dimensional terms and use formulas to convert scanplane pixels (picture elements) into 

area units. For the scan cones substantially equivalent to the 3D conic array 340, the 

algorithms are expressed in three-dimensional terms and use formulas to convert voxels 
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(volume elements) into volume units. The heat-filtering algorithm is determined one of the 
simplest PDE for image processing is the heat equation as shown in equation El : 



du _d u d u 
dt~ dx^ /dy^ 



E 1 



where u is the image being processed whose initial value is the input image: - 1 . 

The above PDE can be discretized using finite differences to give the following 
iterative scheme of equation E2: 



'J I. J 



{u 1 — 2m' • -\- u , ) I — 2u + w' ,) 

V**; + l,y , '^"1,7. ^ ',7 + 1 I, J ■ 



E2 



10 



which can be simplified as shown in equation E3; 



E3 



where k is the time step, and ^ is the value of the pixel at coordinates (/j) at time 



point (iteration) t. 



15 



This heat equation results in a snioothing of the image substantially equivalent to the 



Gaussian filtering of the image. The larger the number of iterations the more the input image is 
; smoothed or blurred. 

Shock filters are another type of PDE-based filters that were originally proposed by 
Osher and.Rudin (Osher and Rudin 1990) for sharpening of signals and images; A one- 
20 dimensional shock filter is given as equation E4; 
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du . ,d u. 



du 



dx 



E4 



where u is the image being processed, — z- is the second partial derivative of u, 

• .' . • ■ . dx ■■ . ■ . 



du 



dx 



': : is the absolute value of the first derivative of u along the x-axis. Here, the signal is modified 
based on the sign of its second derivative and the magnitude of its first derivative. This filter 
5 cannot be discretized using the simple centered finite-difference scheme, 

The result of applying this filter on noiseless blurred data is shown in Figure 5 A. The 
X-axis shows the index coordinates of the signals while the y-axis shows the magnitude of the 
signals. The solid line is the original noiseless signal which was blurred with a Gaussian filter 
algorithm. The blurred signal is shown as the dotted line. A shock filter was then applied to 
10 the blurred signal and that results in the signal shown in light shaded line - it is virtually the 
same as the original signal and the edge is restored from the blurred signal. 

Applying this filter directly to noisy data,, however, does not work and leads to many 

- • • .- ■ . • ■ • . ■ ' • . • .• ' ■ . ■ ' = ' . ■ ■ ^ . ■ 

possible false edges. The combination of heat and shock filters would works very well at 
reducing noise and sharpening of a signal while preserving the significant discontinuities on 

15 a signal. First, the hoisy signal is filtered using the ID heat filter (Equation El), which results 
in the reduction of noise and smoothing of edges. This step is followed by a shock-filtering 
step (Equation E4), which results in the sharpening of the blurred signal. A ID simulated 
example of this conibina;tion heat/shock filter is shown in Figure 5B. The solid line, is the 
original noiseless signal while the dotted dark line is the noisy blurred signal. The result of 

20 heat filtering the noisy signal by the gray dashed line. Finally, the result of shock filtering 

after heat filtering is shown by the gray dot-dashed line. Note that on the last signal, the 

. " ■ . . - - .■ - - " . . i( ' • ■ . ' . _• 

significant edge has been preserved while the noise has been removed. 
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In two dimensions, the shock filter PDE specified by Osher and Rudin is: 



du 
'dt 



= -i^(£(M))||Vw 



E5 



. where ^(m) is the Laplaciaii of the image (u, F,) is a function of the Laplacian, and 

'-""I ' " * - • ' - ■ 

Vw is the gradient magnitude. For the Laplacian, the second derivative in the direction of the 



gradient is used: 



E6 



- For the function to avoid placing spurious edges at points with small gradient 
values as would occur by choosing the sign by selecting the sign of the Laplacian in the 1 -D 
case, the following function is used instead: : / ; 



10 



Fii(uy) = 1 , if f («) > 0 and IV w 



. =-l, if £(w) <0 and Vw 



■ ;\. V. ^ .y^ ' .— 0, otherwise ;- . '.. .•■ 

where Y is a tlu-eshold on gradient values. ^ : . ■ 

Once again, as in the' ID case, a combination of the heat and shock filters works very 
15 well at reducing noise and sharpening. of a signal while preserving the significant edges on an 
image. Figure 2-image 208 shows the; results of applying the heat and shock filter on real 
ultrasound data. . 

The point-positioning algorithm uses two different methods for locating the initial edge 
points. One method looks for edges along radial lines and the other looks for edges along At 
20 mode lines. The first rnethod requires the user to specify a center point from which radial lines 
are drawn for an image, and the second rnethod is completely automatic. 
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In the first method, to locate the initial edge points on an image, the image is radially sampled 
from the user marked center location. The radial lines are drawn at substantially in 30-degree 
increments. On each radial line, edge points are found by using the values obtained from the 
spatial gradient algorithni. A first point on each radial line whose gradient magnitude exceeds a 
5 pre-set threshold and which represents a transition from a dark area to a bright area. If no such 
point is found, that radial line is skipped. As shown in Figure 2 image 210, the points are 
assigned to the edge locations. 

: In the second method, the A-mode lines starting fi'om the transducer position are moved 

outwards. One-dimensional gradients and pixel intensity values / are used along these A-mode 

, .. ■ ■ .■ " . I- - - ' . - - ' , ' ^ . ■ 

10 lines to determine fluid wall locations. A hysteresis threshold is appUed to the . gradient values 

to detenriine the leading wall and the trailing wall of the structure. 

Given two points on an image, pi and p2, an optimal path finding algorithm finds the 

"best" path between the two points where "best" is the one that has the least cost. The cost 

fijnction for a path can be defined to depend upon many factors such as, how close the path 
15 passes through edges on an image, the direction of the path, how smooth the path is, and what 

the length of the path is. A dynamic programming algorithm can be then used to locate the least 

cost path between pi and p2. 

The total cost fiinction to be a weighted is deifined by a sum of 3 terms: 



dp) = aC, (p)^ PC, (p) + }CXp) 



E7 



20 



where d^p^andy are niimerical values and the three cost terms are: 



1. Edge Distance Cost C^(p): This cost is defined such that paths closer to 
high gradient points, i.e., edges, are favored over other paths. 
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where // is the pixel intensity at point i, V/. is the gradient magnitude 



of point i, and is a small constant. Since this term is the sum of 
inverse image gradients for each pixel on the path - this cost term 

■ ■ ' ■ • ■- - ■ ■■ ' " ■ ' ' . . ■'■ . T-r ' ' ■ '■ 

implicitly favors shorter paths over longer paths. 



2. Path Direction Cost: Since the path between p and p isa directed graph, a 

cost can be defined based on the orientation of the graph. The cost is defined to 

■. ' - ' - - ' • ' ' ■ ' ■ - " ■ ■ ■ ... .... 

favor paths for which pixels on the inside of the path are dark and the pixels on 
the outside are bright (See Figure 6). The PDE formula for calculating the path 
direction cost is provided in equation E9. 



where 7^'" and 7, are the image intensities on the inside and the outside of the . 

path connecting point / and L\ as shown in Figure 6. Figure 6 depicts, at the 
pixel level, how the optimal path algoritimi is applied. Direction of path 
proceeds from point 610 (pi.{) to point 614 (pi) between the outside pixel 602 

) and the inside pixel 606 {!'" ). Given two points on the image, 610 

and Pi 614, an optirnal path finding algorithm finds the "besf path between the 
two points where "besf ' is the one that has the least cost and bisects the <; Thus 



^ I 
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pixel 606 is the image intensity of the inside and dark pixel and /f"' 610 is 
the image intensity outside and bright pixel along a path bisecting the 606 
and /f"' 602 pixels and connecting pixels points / 614 and L\ 610 adjacent to 
each side of the 606 pixel. 

3. Previous Contour Distance Cost: When trying to find paths in 3D images, a 
constraint upon the paths on the current, image are imposed and are close to a 
contour on a previous image. Therefore, a definition of a cost based on the 
distance of a point on the current path to a contour from a previous image can 
be made as according to equation E 10: 



p2 



E 10 



where P is the previous contour and Df is the distance of point i from 

the closest point on the previous contour. This distance can be easily evaluated 
\ ': using the Euclidean distance transform fionction computed for the previous 
■ contour. "■ \ - ' \: \ 
. . A dynamic programming algorithm was used to find the least cost path between the 
initial edge points. The two steps typically involved in this are: (1) Computing the cost of the 
path from the starting point to each point on the image - this cost image (also known as the 
minimum action surface) using the fast-marching algorithms; (2) Locating the path coimecting 
the end point to the starting point by gradient descent on the minimum action surface. 

This multi-step segmentation algorithm of the present invention serves to establish a 
determination of the best path between two points and makes the segmentation much more 
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robust. When the t>yo points are quite far apart on the image, it still locates a good path that 
passes through the significant edges and completing, the areas with n^ 

Figure 2 image 212 shows the optimal path found using this algorithm connecting the 



initial edge points. 



s 



The first image processed was the central image that had the fewest missing lateral 
boundaries. No previous contour constraint was used for tWs central image. 

From this central image, die process was moved up to process all the images after the 
central image through the end of the sequence and then moved down to process all the images 
before the cenh*al image through the beginning of the isequence. 



10 



After the contours on all the images have been delineated, the volume of the segmented 



structure is computed/ The contour inside the polygon is filled using a polygon filling 
- algorithm to compute the volume. Then, the intensity value of a central pixel is determined 
using trilinear interpolation (where intensity values from six three-dimerisionail neighbors are 
used) to create a Cartesian volume fi"om the polar coordinate images. Pixels on the boundary 
15 are assigned fractional values depending on the partial volume coverage. Finally, the sum of all 
the pixels in the reconstructed volume is multiplied by the voxel resolution to get the volume. 
The segmentation algorithm was made on images firom a bladder phantom and on 

^.images of the bladder firom several subjects. Boundaries were delineated on all of the images in 

/ • ■ " ■ - - ■ ■ . . -' - ■ , - ■ ■ - " ' ■ -■ ' ■ • . 

* the 3D dataset using the same center point for all images in the 3D sequ^^^ 



20 



„ The segmentation algorithm was first tested on a bladder phantom with a known 



volume. From, the bladder phantom 24 plaries of data are acquired. The boundaries detected on 
this phantomi are shown in Figure 7. 
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While the true volume of the bladder phantom is known to be 495 mL, the volume 
computed using our algorithm is 483 mL, or 97.6% accurate. 

The accuracy for clinical measurements fared well. The segmentation algorithms as 
applied to 12 scan planes of real bladder images from five different subjects (3 females and 2 
5 males) are shown in Table 1. The volume of urine was measured using a UroFlow device. The 
segmentation algorithm of the present invention was applied to images to measure volumes and 
compared to the urine volumes measured by the UroFlow device. The volume measured by the 
Urdflpw device added to the post-Uroflow volume is considered the true bladder as a standard 
to compare to the volume determined by the segmentation algorithms applied to the 12 bladder 



1 0 scan images from each patient. 



The table below shows the raw results on the five subjects- 
Table 1 : Bladder volume computational 



results 



Subject 


True 


Segmentation 


Accuracy 


ID 


Bladder 


Algorithm (mL) 


r 

.. .(%) , 




Volume 




1 t 




:(mL) 






102 


391 


■ . ■ ■ 367 . 


93.9% 


130 


238 


- 232 • 


■ . . 97.5% 


138 


651 

* 


642 


98.6% 


,148 


469 


; .434 


92.5% 


162 


169 


104 


. : 61.5% 
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.Accuracy is defined as a % ratio and is calculated by taking the volume determined 
from the Segmentation Algorithm applied to the 12 image scans and dividing it by the True 
Bladder Volume, then multiplied by 100. The accuracy range of the 12 scan plane 
measurernents to measured urine volumes is 6 J. 5% to 98.6%. The median is 93.9%. The mean 
accuracy for the five measureinents is 88.8%. 

Figure 8 and Figure 9 respectively show the delineated bladder boundaries of 12 scan 
plane images for patients 138 and 148 subjected to the segmentation algorithm pf the present 

» - -- ^ _ ",■ 

invention. The segmentation algorithm set boundaries very reasonable boundaries except at a 
few planes irt that the initial edge point is outside the expected location. 



10 



\The segmentation algorithm of the present invention involves fitting minimum cost 



contours to known boundary points, The cost function is defined very similarly to snakes and 



V ; other deformable models; however, instead of modifying an initial contour, our algorithm finds 
V ah optiriial path between two known points on the boundary. The advantages compared to 
deformiable models are: (a) The initial set pf points needs to be very small and you do not need 
15 an entire initial contour, (b) The algorithm is more robust with respect to missing boundaries 
since missing boundaries are simply handled through an least cost path through them, (c) The 
algorithm is not iterative and alvyays finds a unique global minimum for the cost fiuiction. 

While the preferred embodiment of the invention has been illustrated and described, as 

, . - - -I ■ - - ' . - . - • - A - * - . ^ . - ■ . 

■ noted above, many changes . can be made without departing from the spirit and scope of the 
20 invention. For example, other uses of the invention include determining the areas and volurnes 
' . of the prostate, heart, aiid other organs and body regions of clinical interest. Accordingly, the 
scopeof the invention is not limited by the disclosure of the preferred embodiments. ' 
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